Evaluating conserved domains and motifs of decapod gonadotropin-releasing hormone G protein-coupled receptor superfamily.

G protein-coupled receptors (GPCRs) are an ancient family of signal transducers that are both abundant and consequential in metazoan endocrinology. The evolutionary history and function of the GPCRs of the decapod superfamilies of gonadotropin-releasing hormone (GnRH) are yet to be fully elucidated. As part of which, the use of traditional phylogenetics and the recycling of a diminutive set of mis-annotated databases has proven insufficient. To address this, we have collated and revised eight existing and three novel GPCR repertoires for GnRH of decapod species. We developed a novel bioinformatic workflow that included clustering analysis to capture likely GnRH receptor-like proteins, followed by phylogenetic analysis of the seven transmembrane-spanning domains. A high degree of conservation of the sequences and topology of the domains and motifs allowed the identification of species-specific variation (up to ~70%, especially in the extracellular loops) that is thought to be influential to ligand-binding and function. Given the key functional role of the DRY motif across GPCRs, the classification of receptors based on the variation of this motif can be universally applied to resolve cryptic GPCR families, as was achieved in this work. Our results contribute to the resolution of the evolutionary history of invertebrate GnRH receptors and inform the design of bioassays in their deorphanization and functional annotation.

Application of the random forest algorithm to Streptococcus pyogenes response regulator allele variation: from machine learning to evolutionary models.

Group A Streptococcus (GAS) is a globally significant bacterial pathogen. The GAS genotyping gold standard characterises the nucleotide variation of emm, which encodes a surface-exposed protein that is recombinogenic and under immune-based selection pressure. Within a supervised learning methodology, we tested three random forest (RF) algorithms (Guided, Ordinary, and Regularized) and 53 GAS response regulator (RR) allele types to infer six genomic traits (emm-type, emm-subtype, tissue and country of sample, clinical outcomes, and isolate invasiveness). The Guided, Ordinary, and Regularized RF classifiers inferred the emm-type with accuracies of 96.7%, 95.7%, and 95.2%, using ten, three, and four RR alleles in the feature set, respectively. Notably, we inferred the emm-type with 93.7% accuracy using only mga2 and lrp. We demonstrated a utility for inferring emm-subtype (89.9%), country (88.6%), invasiveness (84.7%), but not clinical (56.9%), or tissue (56.4%), which is consistent with the complexity of GAS pathophysiology. We identified a novel cell wall-spanning domain (SF5), and proposed evolutionary pathways depicting the 'contrariwise' and 'likewise' chimeric deletion-fusion of emm and enn. We identified an intermediate strain, which provides evidence of the time-dependent excision of mga regulon genes. Overall, our workflow advances the understanding of the GAS mga regulon and its plasticity.

Lessons Learnt From Using the Machine Learning Random Forest Algorithm to Predict Virulence in Streptococcus pyogenes.

Group A Streptococcus is a globally significant human pathogen. The extensive variability of the GAS genome, virulence phenotypes and clinical outcomes, render it an excellent candidate for the application of genotype-phenotype association studies in the era of whole-genome sequencing. We have catalogued the distribution and diversity of the transcription regulators of GAS, and employed phylogenetics, concordance metrics and machine learning (ML) to test for associations. In this review, we communicate the lessons learnt in the context of the recent bacteria genotype-phenotype association studies of others that have utilised both genome-wide association studies (GWAS) and ML. We envisage a promising future for the application GWAS in bacteria genotype-phenotype association studies and foresee the increasing use of ML. However, progress in this field is hindered by several outstanding bottlenecks. These include the shortcomings that are observed when GWAS techniques that have been fine-tuned on human genomes, are applied to bacterial genomes. Furthermore, there is a deficit of easy-to-use end-to-end workflows, and a lag in the collection of detailed phenotype and clinical genomic metadata. We propose a novel quality control protocol for the collection of high-quality GAS virulence phenotype coupled to clinical outcome data. Finally, we incorporate this protocol into a workflow for testing genotype-phenotype associations using ML and 'linked' patient-microbe genome sets that better represent the infection event.

In silico characterisation of stand-alone response regulators of Streptococcus pyogenes.

Bacterial "stand-alone" response regulators (RRs) are pivotal to the control of gene transcription in response to changing cytosolic and extracellular microenvironments during infection. The genome of group A Streptococcus (GAS) encodes more than 30 stand-alone RRs that orchestrate the expression of virulence factors involved in infecting multiple tissues, so causing an array of potentially lethal human diseases. Here, we analysed the molecular epidemiology and biological associations in the coding sequences (CDSs) and upstream intergenic regions (IGRs) of 35 stand-alone RRs from a collection of global GAS genomes. Of the 944 genomes analysed, 97% encoded 32 or more of the 35 tested RRs. The length of RR CDSs ranged from 297 to 1587 nucleotides with an average nucleotide diversity (π) of 0.012, while the IGRs ranged from 51 to 666 nucleotides with average π of 0.017. We present new evidence of recombination in multiple RRs including mga, leading to mga-2 switching, emm-switching and emm-like gene chimerization, and the first instance of an isolate that encodes both mga-1 and mga-2. Recombination was also evident in rofA/nra and msmR loci with 15 emm-types represented in multiple FCT (fibronectin-binding, collagen-binding, T-antigen)-types, including novel emm-type/FCT-type pairings. Strong associations were observed between concatenated RR allele types, and emm-type, MLST-type, core genome phylogroup, and country of sampling. No strong associations were observed between individual loci and disease outcome. We propose that 11 RRs may form part of future refinement of GAS typing systems that reflect core genome evolutionary associations. This subgenomic analysis revealed allelic traits that were informative to the biological function, GAS strain definition, and regional outbreak detection.

Phylogenomic history of enigmatic pygmy perches: implications for biogeography, taxonomy and conservation.

Pygmy perches (Percichthyidae) are a group of poorly dispersing freshwater fishes that have a puzzling biogeographic disjunction across southern Australia. Current understanding of pygmy perch phylogenetic relationships suggests past east-west migrations across a vast expanse of now arid habitat in central southern Australia, a region lacking contemporary rivers. Pygmy perches also represent a threatened group with confusing taxonomy and potentially cryptic species diversity. Here, we present the first study of the evolutionary history of pygmy perches based on genome-wide information. Data from 13 991 ddRAD loci and a concatenated sequence of 1 075 734 bp were generated for all currently described and potentially cryptic species. Phylogenetic relationships, biogeographic history and cryptic diversification were inferred using a framework that combines phylogenomics, species delimitation and estimation of divergence times. The genome-wide phylogeny clarified the biogeographic history of pygmy perches, demonstrating multiple east-west events of divergence within the group across the Australian continent. These results also resolved discordance between nuclear and mitochondrial data from a previous study. In addition, we propose three cryptic species within a southwestern species complex. The finding of potentially new species demonstrates that pygmy perches may be even more susceptible to ecological and demographic threats than previously thought. Our results have substantial implications for improving conservation legislation of pygmy perch lineages, especially in southwestern Western Australia.

In silico characterisation of the two-component system regulators of Streptococcus pyogenes.

Bacteria respond to environmental changes through the co-ordinated regulation of gene expression, often mediated by two-component regulatory systems (TCS). Group A Streptococcus (GAS), a bacterium which infects multiple human body sites and causes multiple diseases, possesses up to 14 TCS. In this study we examined genetic variation in the coding sequences and non-coding DNA upstream of these TCS as a method for evaluating relationships between different GAS emm-types, and potential associations with GAS disease. Twelve of the 14 TCS were present in 90% of the genomes examined. The length of the intergenic regions (IGRs) upstream of TCS coding regions varied from 39 to 345 nucleotides, with an average nucleotide diversity of 0.0064. Overall, IGR allelic variation was generally conserved with an emm-type. Subsequent phylogenetic analysis of concatenated sequences based on all TCS IGR sequences grouped genomes of the same emm-type together. However grouping with emm-pattern and emm-cluster-types was much weaker, suggesting epidemiological and functional properties associated with the latter are not due to evolutionary relatedness of emm-types. All emm5, emm6 and most of the emm18 genomes, all historically considered rheumatogenic emm-types clustered together, suggesting a shared evolutionary history. However emm1, emm3 and several emm18 genomes did not cluster within this group. These latter emm18 isolates were epidemiologically distinct from other emm18 genomes in study, providing evidence for local variation. emm-types associated with invasive disease or nephritogenicity also did not cluster together. Considering the TCS coding sequences (cds), correlation with emm-type was weaker than for the IGRs, and no strong correlation with disease was observed. Deletion of the malate transporter, maeP, was identified that serves as a putative marker for the emm89.0 subtype, which has been implicated in invasive outbreaks. A recombination-related, subclade-forming DNA motif was identified in the putative receiver domain of the Spy1556 response regulator that correlated with throat-associated emm-pattern-type A-C strains.

In silico prediction of the G-protein coupled receptors expressed during the metamorphic molt of Sagmariasus verreauxi (Crustacea: Decapoda) by mining transcriptomic data: RNA-seq to repertoire.

Against a backdrop of food insecurity, the farming of decapod crustaceans is a rapidly expanding and globally significant source of food protein. Sagmariasus verreauxi spiny lobster, the subject of this study, are decapods of underdeveloped aquaculture potential. Crustacean neuropeptide G-protein coupled receptors (GPCRs) mediate endocrine pathways that are integral to animal fecundity, growth and survival. The potential use of novel biotechnologies to enhance GPCR-mediated physiology may assist in improving the health and productivity of farmed decapod populations. This study catalogues the GPCRs expressed in the early developmental stages, as well as adult tissues, with a view to illuminating key neuropeptide receptors. De novo assembled contiguous sequences generated from transcriptomic reads of metamorphic and post metamorphic S. verreauxi were filtered for seven transmembrane domains, and used as a reference for iterative re-mapping. Subsequent putative GPCR open reading frames (ORFs) were BLAST annotated, categorised, and compared to published orthologues based on phylogenetic analysis. A total of 85 GPCRs were digitally predicted, that represented each of the four arthropod subfamilies. They generally displayed low-level and non-differential metamorphic expression with few exceptions that we examined using RT-PCR and qPCR. Two putative CHH-like neuropeptide receptors were annotated. Three dimensional structural modelling suggests that these receptors exhibit a conserved extracellular ligand binding pocket, providing support to the notion that these receptors co-evolved with their ligands across Decapoda. This perhaps narrows the search for means to increase productivity of farmed decapod populations.